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A mimetic discretization of the Abelian Chern-Simons theory is presented. The study relies on 
the formulation of a theory of differential forms in the lattice, including a consistent definition of the 
Hodge duality operation. Explicit expressions for the Gauss Linking Number in the lattice, which 
- - - correspond to their continuum counterparts are given. A discussion of the discretization of metric 

structures in the space of transverse vector densities is presented. The study of these metrics could 
. serve to obtain explicit formulae for knot an link invariants in the lattice. 
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I. INTRODUCTION 



The relationship between quantum Chern-Simons theory and knot theory is an active field of research both in physics 
and mathematics P, The key contribution in this field, due to E. Witten, stated that the vacuum expectation 
value of the Wilson loop in the Chern-Simons theory yields knot invariants related with the Jones polynomial 1, 3]. 

Besides their interest from the mathematical side, knot invariants play an important role in the loop representation 
of quantum gravity [^,0] based in Ashtekar variables jdj . In this formulation, wave functionals are just knot- invariants, 
hence, it is useful to have analytical expressions for knot-invariants that can be properly handled through the machinery 
of loop-calculus. A systematic way to obtain expressions of this kind consists on considering the perturbative expansion 
of the Wilson loop in the quantum Chern-Simons theory 0, Under very general conditions, it can be shown that 
each order in this perturbative expansion is a knot invariant. Hence, following this method one can get, in principle, 
an infinity of such invariants. It should be said that there are also "knot" and "link" invariants associated to the 
spin- networks formulation of quantum gravity Q, which is the modern setting in which loop quantum gravity is 
studied. 

There also exists another method to obtain knot-invariants, which deals with the classical field formulation. This 
f — ■ method consists on considering the on-shell action of a diffeomorphism-invariant theory coupled to external sources 
I with support on a curve (or several curves, if one desires to consider link- invariants instead of knot-invariants). The 
(<~^ . metric independence of the theory guarantees the topological character of the on-shell action, which yields analytical 
^-H ' expressions for knots (or link) invariants. This method was first considered for abelian theories ilOi|, and then extended 
CS| . for non-abelian ones 

■ Regardless one uses the quantum or the classical method, the link- invariants (i.e., those invariants involving several 
^ I loops) that rise from the Chern-Simons theory are well defined (provided that the curves involved do not intersect 
each other). This is not always the case for knot- invariants, which can suffer from ambiguities (ref), as the Gauss 



self-linking number exemplifies [12| . A way to cure these ambiguities consists on choosing a " framing" : instead of 
considering a curve, one takes a strip whose width will be taken to vanish after the entanglement properties are 
calculated. Another way of handling the ambiguities is through the discretization of the knot or link invariants, since 
the lattice can provide a natural framing. On the other hand, a discretization of topological field theories or loop 
quantum gravity can be also useful for calculational purposes. 

There are many ways to discretize a physical theory, since there are many forms to convert differential equations 
into finite difference equations, or to set a sensible discretization of vectorial operators such as curl, gradient and 
divergence. The line, surface and volume integral of vectorial fields admit different discretizations as well. When 
the discretization is made in such a way that most of the significative mathematical properties (in the continuum) 
of the objects involved are preserved in the lattice, one says that it is "mimetic". There is a vast literature [Ijl 
about mimetic discretizations in the framework of classical mechanics and electromagnetism. Also, the Hamiltonian 
approach of QCD in the lattice [U] implicitly deals with a mimetic discretization of the Yang-Mills theory, although 
this fact is seldom stressed. 

Some years ago, a "consistent discretization" of the linearized theory of General Relativity, both in the traditional 
metric variables and in the Ashtekar variables was considered, yielding a mimetic discretization of the theory [15j. A 
mimetic discretization of the Abelian Chern-Simons theory in a lattice was proposed by Frhlich and Marchetti jll] 
and studied by several authors 17]. A mayor problem related to the discretization of theories with differential forms 
consists on the mimetization of the Hodge dual operation. In reference [l^ another mimetic discretization of the 
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Chern-Simons theory was considered, and using the quantum method explained above the authors obtained a formula 
for the self-linking number in the lattice. There was argued that the framing problem is related to the discretization 
of the Hodge dual operation. Hence, a prescription for this operation was given in that reference, but it is knot 
completely mimetic, since taking the dual twice does not lead back to the starting form. 

In this article we study a mimetic discretization of the Abelian Chern-Simons theory and discuss how to obtain 
knot invariants in the lattice from the evaluation of the discrete on-shell action. This is achieved by means of the 
introduction of a formalism for the mimetic discretization of differential forms in the three-dimensional cubic lattice, 
which includes a natural formulation of Hodge duality. We write down the Chern-Simons action coupled to external 
currents associated to curves in the lattice (which are also mimetically discretized) , obtain the equations of motion, 
solve them, and calculate the action on shell, which produces formulae for knot (or link) invariants in the lattice. 
Also, we address (to some extent) the problem of taking the analytical expressions for some higher order knot and 
link invariants already known in the continuum, and producing discrete formulae for them in a mimetic form. To 
this end, we elaborate further on the study of the "metrics" in the space of transverse functions studied in reference 
p^ . These metrics can be seen as building blocks in the construction of invariants, both in the continuum and in the 
lattice. 



II. DIFFERENTIAL FORMS IN THE CUBIC LATTICE 

We shall denote hy v, I and p, the sites, oriented links and oriented plaquettes of the cubic lattice, respectively, . I 
(p) denotes the link (plaquette) that coincides with / (p) but has opposite sense. The elementary cubes, plaquettes, 
links and sites are the the g-cells of the lattice, with g = 3, 2, 1, respectivelly. L means the distance between nearest 
neighboring sites. Sometimes it will be convenient to label the g-cells using sites and unitary vectors as follows. Given 
a right handed orthonormal basis {ei, 62, 63} formed by vectors parallel to the links of the lattice, we can write, for 
instance, v' = v + Le, meaning that the sites v' and v are nearest neighbors separated by a link in the direction of e. 
We shall also use the notation 

Li = Lei and L123 = L(ei -f 62 -I- 63) . (1) 

If the link I starts at v and ends at w' = w -I- L e, we shall also set I — L{v, e). With this notation, the opposite link 
is Z = L{v' , — e). The site where I starts will be denoted as .1 and the final one will be I. We say that a link I pointing 
along one of the basis vectors is a positive one, and we write Z > 0. Consider a plaquette coplanar with and Sj, 
whose orientation is given by x ij. Let v the site of p having lower values of its coordinates. Then we say that 
p = P(u, ei,ej). A plaquette whose normal vector points along some of the e.^ is a positive one, and we write p > 0. 

A discrete 0-form is a function taking values on the sites of the lattice. If A is a 0-form, A„ is the value it takes at 
site V. The domain of a discrete 1-form, instead, is the set of links of the lattice. Let A be a 1-form and I = L{v, ii) a 
link; the evaluation of ^ at Z will be indistinctly denoted by Ai = Ai{v). We shall also say that Ai{v) is the component 
along ii of the vector field A evaluated at v. We ask the 1-forms to obey 

At=-Ai, (2) 

which is the lattice version of the "continuum" equation A ■ (— e^) = —A ■ (e^). A discrete 2-form takes values on 
the plaquettes: Bp — Bij{v) means the evaluation of the 2-form B on the plaquette p = P{v^ e^, Cj). We impose the 
condition 

Bp = —Bp or equivalently Bji{v) — —Bij{v) . (3) 

Also, we set Bii{v) = 0. As in the continuum, discrete 2-forms are antisymmetric. Finally, we define a 3- form as a 
function taking values at the cubes c of the lattice. If F is a 3-form, we write 

Tc = Fi 2 3 (w) = Tj J fe (w) = -Tkj iiv) (4) 

Here, v is the site with lower coordinate values belonging to c, and fc} is a cyclic permutation of {1,2,3}. We 
also set Tijk{v) = whenever any pair of indexes take the same value. 

Given a 0-form f{v), we define its discrete partial derivative in the direction of as 



d,f{v)=f{v + U)-f{v). 



(5) 
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It is easy to see that 

djdifiv) = didjfiv) , (6) 

and that when the lattice spacing tends to zero, the discrete derivative tends to the continuum derivative. 
We define a g-form ^ as a function over the g-cells. Its exterior derivative {d(j)) is a (g' + l)-form defined by 

(#)(9+l)-cell = X/ ^g-ce\h (7) 

g-cellea(g+l)-cell 

where d{q + l)-cell is the boundary of the {q + l)-cell. Let us see how this formula looks in detail. For a 0-form A we 
have 

{dX)i = \i. - A., . (8) 

If / = L{v, Ba) we have 

{dX)L(v,e^) = X{V + L„) - \{V) = daX{v) . (9) 

Now, let A be a 1-form. Its exterior derivative is given by 

{dA)p = J2Ah (10) 

where the links I belonging to the edge of p are taken accordingly to the right hand rule with respect to the orientation 
of p. U p = P{v, e,, Cj), then 

{dA)ij{v) = Ai{v) - Ai{v + hj) + Aj{v + U) - Aj{v) , 

hence 

idA),,iv) = d,Aj{v) - djA,iv) . (11) 

Next, take 5 to be a 2-form. We have 

{dB), = Y,Bp (12) 

where the plaquettes p that bound the cube c are oriented outwards. If v is the site of c having lesser values of their 
components, then 

(dB)c = -B3l{v)+B3l{v+L2)-Bi2{v) + Bi2{v+I.3)-B23{v)+B23{v+U) = d2B3l{v) + d3Bi2{v) + diB23{v) . (13) 

Therefore we have 

{dB)ijk{v) = d(iB^k)o{v) = ddBefiv) , (14) 

where dfj^ stands for the antisymmetrized product of three Kronecker's deltas. Finally, the exterior derivative of a 
3-form is taken to be zero. Summarizing, if ^ is a g-form, its exterior derivative obeys 



1 

It can be easily shown that 

d# = (16) 

for any discrete form (j). Thus, we have a mimetic version of differential forms and exterior derivative in the cubic 
lattice. 

There exist many ways for defining the Hodge dual #a; of a form oj in the lattice. We shall define it by relating 

forms in the "direct lattice" R with forms in the "dual lattice" K'', which is defined, as usual, in such a way that its 
sites coincide with the centers of the cubes of the former one R. Hence, links of the dual lattice go across plaquettes 
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of the "direct" lattice, while plaquettes of the dual lattice are crossed by links of the "direct" one. Our total space 
will be the union of both K and M.'^. If we want to stress that an elementary cell belongs to a lattice which is the dual 
one of another one, we shall attach the superscript *. Given the 0-form A, we define its Hodge dual #A as the 3-form 
given by 

(#A)e. = A, (17) 

where c* is the cube whose center coincides with the site v (c* belongs to the "dual lattice" if v belongs to the "direct" 
one). The components {4f^)ijk{v*) satisfy 

{#X)ijk{v*) = eijkX{v) with u* = u - ^ L123 • (18) 

For a 1-form A, we define its Hodge dual #A as the 2-form that satisfies 

(#A)p. = Ai (19) 
where p* is the plaquette dual to the link /. It can be seen that 

{#A)ij{v*)=eijkAk{v) with = w + ^ L123 - Lj - L^- . (20) 

Now let B be a 2-form. Its Hodge dual #5 is defined as the 1-form satisfying 

i#B)i, = Bp (21) 
where I* is the link dual to the plaquette p. We have 

{#B)i{v*) = ^eijkBjk{v) with V* =v + ^U23-'Li- (22) 
Finally, given a 3-form F we define its Hodge dual #F to be the 0-form which fulfils 

(#r)„. = r, (23) 

with V* being the site dual to the cube c. One also has 

{#r){v*) = ^Sijknv)ijk with v*=v+^U23. (24) 
These definitions and results may be summarized as follows. If is a g-form, its Hodge dual satisfies 

(#0)ii-i3-,(«*) = ^_^ii-i3-,3i-j<,(t^h-jA^) (25a) 

where 

^* = + I (ei + 62 + 63 - 2e,, 2e,3_ J . (25b) 

It should be noticed that v* belongs to the lattice dual to that where v lives. An immediate consequence of the 
definition of Hodge duality here presented is that for any discrete form one has 

##0 = </- (26) 

just as it happens in the continuum framework. It should be underlined that there exist discretizations of the Hodge 
duality operation that do not "mimic" this property of the continuum. 

We also introduce the coderivative of a discrete g-form, which we define as the {q — l)-form given by 

{6^) = {-iy#d#(t>. (27) 

A careful application of the above definition leads to the result 

((5(/))(q_i).cell = X] (Pq-cell (28) 

Q-ccU/ 
(q-l)-celleOq-cell 
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where dq-ce\\ is the boundary of the g-cell. To obtain the components of the coderivative, we find it convenient to 
introduce the "overhned" partial derivative 

d^fiv) EE f{v) - f{v - ha) - dafiv - LJ. (29) 

Then it can be shown that 

= -9fc0fcii...i,_i(u -Lfe) = -9fe0fcji...i,_i(w) . (30) 
Successive apphcations of the coderivative yield 

5S(j) = 0. (31) 

Next, we show how to discretize integration by parts in a mimetic form. In the formula that follows, the summation 
is over all the sites of the infinite cubic lattice. We have 

At, dafJ'v = At,(/it,+L„ — ^J'v) = ^^(At,-L„ — K})^J'V — — ^^('^u ^ At,_L„)/it, , 

V V V V 

hence 

^^XvdaHv = -^^{daXv-I^Jfiv = - ^(9aA„)/it, . (32) 

V V V 

It will also be useful the following definition. Given two (/-forms in the lattice, their inner product will be 

(0, V') = '0)r + V')R<i (33) 

where 

((/), ■0)h = ^ '/'g-coU'0g-CGll = (V", 0)m (34) 

ij-coll>0/ 
9-C0II6R 

(0, V')Rd= ^ (f>q-ccll-lpq-ccl\ ^ {-lp,(p)Rd. (35) 
(j-ccll>0/ 

(j-coiieR'* 

Here, g-cell > means that the cells are taken with just one of the two possible orientations. In components, the 
inner product is written as 



and 

('^''^') = ^E<^H ■■■.,(«)<. ..,,(«)■ (37) 

The inner product satisfies the following properties 

(d0,V)R = (</',5V')r (38) 
((/.,#7/;)r = (#0,V)R<i. (39) 

that mimic well known properties of the continuum. From them it is immediate that 

= (</',<5V') (40) 
^). (41) 

To prove psp let us take a g-form and a (g + l)-form ^. Then 



dGR uGR 

-^^afc<?!>»i-»,(w)V'fe*i ■■■»,("). 
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But 



therefore 



In turn, to prove ([39]) consider a g-form and a (3 — g)-form ij:. We have 



(3-g) 
1 1 



with 

w* = ti + ^(ei + 62 + 63 - 2eii 2ei3_ J = w - ^(e^ + 63 + 63 - 26^^ 2ejJ 



On the other hand 



(3-g)! ''■■■''-'^'■■■^'"'^■■■'^-'^ ' (3-g) 
hence we finahy have 

To conclude this section we define the discrete versions of the most common differential operators of vectorial calculus: 
gradient, curl, divergence and Laplacian. Given a scalar field or 0-form A, its discrete gradient Q\ is defined as 

g\ = d\ hence {g\\{v) ^ d,\{v) . (42) 

We define the discrete curl TZA and divergence VA of a vector field or 1-form A as follows. 



7^A = #dA = <5#A, hence {TlAUv*) = -e^Jk{dA)Jk{v) = e^Jk^,Ak{v) (43) 



1 

VA = ~5A, therefore (PA)„ = d.,A,{v - L,) = d.,A,{v) (44) 



where v* ^ v + L123/2 — L^. Since dd = 5 5 = Q one has 

ng ^vn = o, (45) 

which mimics what occurs in the continuum framework. We define the discrete Laplacian of a q-form (j) as the g-form 
given by 

= -(d(5 + Jd)(/), (46) 

and it can be shown that 
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III. ABELIAN CHERN-SIMONS THEORY AND THE GAUSS LINKING NUMBER 

Let us briefly review how the Gauss Linking Number can be obtained through the classical theory in the continuum 
The action of the Abehan Chern-S imons theory, coupled to an external source, can be written down as 

S = J d'x{^e''''Mx)d,Ak{x) - aMx)J'{x)). (48) 

It is metric independent; also, it is invariant under general coordinate transformations provided that the current J^{x) 
transforms as a contravariant vector density of weight 1. The equation of motion is 

e'^'^djAkix) =a.rix). (49) 

For this equation to be consistent, the current must be conserved. This fact, as in Maxwell theory, is deeply related 
to gauge invariance. The equation of motion, being formally Ampere's Law, is solved by the Biot-Savart law 

At: J \x — x'\-^ 

where A is an arbitrary function that carries the gauge freedom mentioned above. If we take the current as the "path 
coordinate" of a closed path (i.e., a loop C) 

X\x)= j> dx''5^{x-x'), (51) 
substitute into the "Biot-Savart" expression and finally in the action, we obtain that the action on-shell yields 

871-/(7 Jc \x-x'\-^ 

which is times the Gauss Linking Number L{C, C), evaluated at C" = C Indeed, the above expression is ill 
defined, since the denominator vanishes whenever x and x' coincide. If we take two different (non-intersecting) curves 
instead of just one, this problem is cured. This is related to the issue of "framing" [l^, which will be discussed latter. 
The Gauss Linking Number can also be written down in terms of line and surface integrals, which will be particularly 
useful for studying the lattice version of this link invariant. To this end, we find useful the following notation. Define 
Greek indexes as "composed indexes" fj, = (ix), comprising a continuous and a discrete one. Also we introduce the 
following summation convention 

W'A^ = Jd'^xB'^A,. (53) 

It is worth noticing that the the path coordinate (which may be defined for open paths as well) obeys 

X(7i -72) =x(7i)+x(72), (54) 
Xil) = -X(7), (55) 

where 71 ■ 72 is the path product and 7 is the path opposite to 7. Also, it verifies the differential constraint 



9ixX"il) / dy'dixd{x - y) ^ ~ I dy'diyS{x - y) = S{x - x.^) ~ d{x - x^.) . 

J y J 'y 

In particular if the path is a loop C one has 

a.x"(C^)=0. (56) 
Hence, there exists a "vector potential" mix{C) of x"(C) obeying 

X"(C) =e^^'=a,,mfc,(C), (57) 
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which is defined up to the gradient of an arbitrary scalar field. To calculate a convenient "vector potential" mix{C), 
take a surface whose boundary is the loop C. From Stokes theorem we have 

0= (fdl-A - [ dS-V X A = (f dy'A.y - [ dSk {v)e^ ^ 'djyAiy 



-I 
-I 



d^xA 



dy'6{x -y)- dSk{y)e''^'djy5{x - y) 



(fdy'S{x -y)- e'^^dj, [ dSk{y)S{x - y) 

Jc JHc 



where Ai^ is any vector field. Hence we can write 

X'-{C) = e'i "dj^Mk.illc) = dj^vi' (Sc) , (58) 
where the particular "vector potential" MkxC^c) and its Hodge dual (or "co-potential") r]'^J{T,c) are given by 

Mfc,(Sc) = / dSk{y)S{x-y) (59) 

V'j{^c)= f e'^HSu{y)6{x-y). (60) 

If u and V parameterize the surface Sc, and its orientation is such that 



dS{y) = -^x —dudv, 



we shall have 



(61) 



e^^''dSk{y)=e'^''skni^^dudv 



Qy[^ yj] 

du dv 



dudv. 



Hence, M(Sc) y vi^c) may also be written as 

Mk,{J:c) = le,jk [ dT^y^5{x-y) 



I Tic 

dT}J 5{x - y) 



with 

In terms of the "vector potential" Mkx{^c) the action on-shell can be cast in the form 

„2 



^onshell 



-— d) dx"^ 
2 Jc 



dSi{y)d{x - y), 



(62) 

(63) 

(64) 

(65) 
(66) 



which provides a simple geometric interpretation of the Gauss linking number: it counts how many times one of the 
curves cuts the surface bounded by the other one (in this case both curves are the same, but we could replace one of 
them by another one in the formula obtained above). 

This discussion can be presented in the formalism of differential forms. Indeed, we shall translate these results to 
the lattice framework taking advantage of the mimetic differential-forms language introduced in the previous section. 
The Chern-Simons action coupled to a current, in the lattice, can be written as 



Scs = ki<A, #dA > +k2 < A,j >=kiY^ Ai#dAi -I- fca ^ Aiji , 

l>0 l>0 



(67) 



9 



where A is a discrete 1-form (the lattice Chern-Simons field) and j is the current. The discrete 3-space comprises 
both the direct and dual lattices, as discussed at the end of the first section. Under lattice gauge transformations 
A + dX, where A is a 0-form, the action becomes 



S' =ki{< A, #dA > +<dX, #dA >) + k2{< A,j > + < dX,j >) (68) 
+ ki < A, 6i^dA > +k2 < A, Sj > . 

Since S^dA = 55^A = 0, we find that the action will be gauge invariant provided dj = 0. This conservation law, in 
turn, implies that j derives from a potential m (or a co-potencial r], such that #m = rj ) 

j = 5=f^m = Srj. (69) 

To study the dynamics of the system, we make arbitrary and infinitesimal variations 6 A of the fields. Then, the action 
varies as 



SS = ki{< A, #d6A > + < #dA >) + k2< 6A,j > . (70) 
The variational principle states that 6S = for any 6 A, hence we have 

ki{6#A + #dA) + k2j = 0, (71) 

and finally 

S#A=-^^, (72) 
To solve this equation, we recall that j derives from a co-potencial r], which yields 

^'-W*-^ (73) 

To proceed further, we miist pick an appropriate choice for the external current j, i.e., one that mimics the loop 
coordinate x"(C) of the continuum. Let us consider a path 7 in the whole discrete space, that is, a succession of 
links belonging to the union of the direct and dual lattices. Then we define the lattice path coordinate xil) as the 
1-form given by 

X/(7) = E(^M'-'5mO- (74) 
I'ej 

Xi (7) counts how many times the path 7 passes along I minus the number of times passing along I. To explore further 
how far xi (7) mimics the former path coordinate let us define, as in the continuum, the contraction of two 1-forms A 
and H by 

H^A^ = HiAi = Hi{v) Ai{v) = {H, A) , (75) 
where the sums are over links and sites of R U K''. The contraction of xil) with the 1-form A is 

l>o I'e-r i>o I'e-r i>o 



hence we finally have 



X^(7)^^ = (X(7),^> = E^''- (76) 



As in the continuum, the path coordinate obeys 



X(7i • 72) = x(7i) + X(72) (77) 
X(7) = -X(7i) (78) 
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Also, a careful calculation shows that the coderivative of xil) is given by 

iSxil))v = ~{T^x{l))v = S.,., - (79) 

hence, for a loop C one has 

Sx{C) - . (80) 

Therefore, x(C) can be written as the coderivative of a 2-forni 77(C) ("co-potential"), or the coderivative of the Hodge 
dual of a 1-form m{C) ("potential") as follows 

X{C) = SviC) - 5#m(C) . (81) 

It should be kept in mind that the space comprises both the direct and the dual lattices. Hence, when we talk about 
a loop we mean a "drawing" that may consist on several closed components, some of them in one of the lattices and 
some more in the other one. 

The potential and copotential associated to xil) can be explicitly constructed as follows. Let be a set of oriented 
plaquettes forming a surface bounded by the loop C. Define the 2-form r]{T,c) as 

In words, ri{^c)p equals 1 (—1) if p G Sc (p G '^c) j and vanishes in any other case. From rjiJlc) we construct the 
1-form 



M(Sc) = #ry(Sc) . (83) 

Notice that Af(Ec) equals ±1 when it is evaluated on the links that "puncture" the plaquettes composing Ec, 
depending on whether these links are oriented as the punctured plaquette or in the opposite sense. In any other case 
M(Ec) = 0. It is not hard to see that 

X{C) = SrjiJ^c) , (84) 

therefore, when one takes x(C) as the external current of the lattice Chern-Simons theory, the solution (|73p of the 
field equation (j72p can be written down directly by replacing 77 by 77(Sc). Then, substituting the solution (1731) into 
the action ([67]) we can finally obtain the following expression for the action on shell 



Sonsheii = - J-(#r7(Sc),x(C)) = -^{Mi^c),xiC)) . (85) 
4fci 4A;i 

This expression defines the Gauss Linking number among two loops Ci, C2 

V{C,,C2) = (M(ScJ,x(C^2)) = ^Af,(ScJXi(C2). (86) 

The right hand side of equation (j86p detects whether or not there are plaquettes belonging to whose dual link (or 
the opposite one) coincides with some link of C2. This amounts to say that ip{Ci,C2) detects if C2 punctures Eci, 
i.e., if C2 links Ci. Notice that this linking is only possible among closed components of Ci and C2 that belong to 
different lattices. If all the components of both loops belong to, say, the direct lattice, the Gauss Linking number will 
always vanish. However, there is a way to detect entanglement among loops in the same lattice, as we shall show next. 
For every site in the direct (dual )lattice, there are eight sites of the dual (direct) one "surrounding" it. Therefore, 
given a drawing in one of the lattices, there are eight different copies of this drawing in the other lattice, which can 
obtained by elementary translations along the oblique eight directions pointing towards the surrounding sites. Then, 
given two loops we define their "smeared" Gauss Linking Number (p'{Ci,C2) as the average of the former linking 
numbers among one of the loops and each copy of the other loop which is obtained by translating it " slightly" in the 
eight possible forms described above. In other words, we set 

1 ^ 

'^'(Ci,C2)--^^(r.Ci,C2), (87) 
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where TaC represents the a'th translation of the loop C (there are 8 such possible translations). It is not difficult to 
prove that 



8 8 

-5](^(r,Ci,C2) = -^^(Ci,r-iC2), 



and that (p'{Ci,C2) is invariant under lattice translations and (discrete) rotations, as one should expect. 

The smeared Gauss Linking Number (f'{Ci,C2) can be seen as a lattice implementation of the idea of "framing" 



in the continuum framework [12|: to avoid ambiguities when studying the self-linking of a loop, one makes a slightly 
translated copy of the loop that do not intersect with it, and calculates the Gauss Linking Number among the original 
and the "framed" loop. Of course, the result depends on the particular choice of framing. In our case, one could say 
that the smeared Gauss Linking Number (p'{Ci,C2) takes into account all the possibilities of framing in the lattice 
on the same footing. An interesting feature of this construction is that ip'{Ci, C2) yields well defined results even in 
the case of intersecting loops. Also, it is interesting to notice that in those cases, unlike the non-intersecting ones, 
the results are not integers. For instance, consider two loops composed by just one closed piece lying in the direct 
lattice. If the loops intersect each other at one site forming a cross, the contribution to (p'{Ci, C2) at that site will be 
i . If the intersection is among a straight segment of one of the loops and a corner of the other one, the result will be 
■|. This is a sensible manner to take into account those limiting cases of linking: one could say, for instance, that a 
crossing is "half a linking. 

IV. METRICS g*^^^ AND g^^jy IN THE CONTINUUM 

The expression for the Gauss Linking Number in the continuum (j52[) suggests the existence of an underlying 
structure that we shall explore in this section: the generalized metrics in the space of transverse functions. This study 
was initiated in the continuum framework in reference [? ] . A vectorial density field is transverse if its divergence 
vanishes. Define the transverse projector as 



where is any function satisfying 



and 



5%^^5fy+d,y<i>l^ (89) 

a„</,;- = 5{x - y) , (90) 
5% = 5{x - y)5,, . (91) 



Since 



d.J^jy^i), (92) 
^Tjy projects a vector field into the space of transverse functions. Also, it acts as the identity operator in that space 

5t^^ T" = if 9„T" = . (93) 

There is a projector for each choice of (j) (later, we shall refer to the relationship among different choices). Finally, we 
have 

hence, S^jy is indeed a projector. We define the covariant g^i, and contravariant g'^'^ metrics in the transverse-forms 
space as 

9ixjy — £ijk4'y : (94) 

g"'y = dtJ{x-y), (95) 
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where 4>t is asked to obey, besides (|90] 



These properties imply that 



hence, the g' s are symmetric 



(96) 
(97) 



9^.0;^ = , (98) 

(99) 

i/^^^'^yJ^", (100) 



and their contraction yields 



The contravariant or inverse metric g^'^ metric obeys 

5„<7"''^ = 9,yff"^^=0, (102) 

and acting on a vector field A„ yields its curl 

g'^^A^ = e'^'^dj^Aux = (7^A)" . (103) 
If is a transverse vectorial density, its potential vector A^ defined by 

T'" = {TZAY or T = 5#A = i^dA , (104) 

can be written down as 

A^^gf^.T". (105) 

In fact 

{TZAr - g'^-'A, = g^^guxT^ = Sl^^T^ = T'^ . 

A choice of 0" determines a choice of metric and a particular solution for the potential vector. Changing 0" amounts 
to making a gauge transformation A' = A + dX. 

Given two transverse vectorial densities Ti and T2 we define their inner product as 

g(ri,r2) =Tf.g^,T2^ (106) 

which is independent of the choice of prescription for 0", due to the transverse character of the T's. The inner 
product (|106[) can also be written in terms of vector potentials. Given — g^'^A^^ and T2 ~ g'^ ^A'p one has 

5(Ti,T2) = Ti^g.^T^ = g^'^A^ g,, g'^^Al = A, „ S^^g'^f' A^ =A\g'^^ A%. 

Taking A = A^—A^ one obtains the Abelian Chern-Simons action. 

Let us study some useful prescriptions. One of them, the "transverse prescription" is given by 



An \x — y\ An \x — y|' 



Since 
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it is immediate to check ([90)) . Also, (|96l) and ([97)1 are fulfilled. The covariant metric corresponding to this prescription 
is given by 



90 ix jy 



A'k\x — 



and it is transverse 



IX ]y 



= 0. 



(109) 



(110) 



Taking the loop coordinates T' s and choosing the metric precisely in this prescription we see that the inner 

product T^g^vT2 just corresponds to the Gauss Linking Number ip[Ci,C2) = X^{'Ci)9fivX'' {C2)- Since the result is 
independent of prescription, it could employed another ones to present different formulae for the Gauss invariant, which 
could be convenient in the lattice framework (the Green's function of the Laplacian that appears in the conventional 
formula in the continuum can not be put in a mimetic fashion). 

The following prescription will be particularly useful in its lattice version. Let hhe a, continuous path going from 
the origin to the spatial infinity. z[t) label the points of the path, r being a parameter and z(0) — the origin. We 
define the curves hy^ obtained by translating h until its initial point coincides with the point y. Then, if z{T,hy) is 
the point of the curve hy corresponding to the parameter value r, one has 



z{t, hy) = z{T)+y. 
It can be shown that the path coordinate x''^ of the curves hy satisfies 

x''=ihiy-.))=x'^''^''Hhy). 
From this construction, let us define the "bundle of paths" prescription as follows 



"h y 



dz^S{x — z) — / dz^5{y — z) 



^^[x'%hy)~x'y{h,)\ 



(111) 



(112) 



(113) 



It can be easily seen that obeys ([96]). On the other hand, using (|112p one has 



y 
1 

2 



hence (^7)) is also fulfilled. Regarding condition (PO)) . one has 



2d, 



■ixVh y 



ii{x+txj) 
'^h y 



dz''dixS{x - z) - dix / dz'5{y-z) = - I dz''dizS{x - z) + diy dz''5{y-z) 

J J hy J 

dz'd,J{x -z)- f dz'd^Jiy -z) = -[-5{x - y)] - [-5{y - x)] - 25{x - y) . 



The covariant metric associated to this prescription is then 



Qhixjy — ^ijk0h 



= '-f[x'''ihy)~x''Hhx)] 



kx 



2 



dz''6{x - z) - / dz''d{y - z) 



(114) 



Regardless which prescription for the metric we employ, it should be stressed that the Gauss Linking number can be 
written as 



Vp(Ci,C2)-x^(Ci)g^,x'^(C2). (115) 

In the "bundle of paths" prescription, the interpretation of ()115|) is as follows. The paths of the bundle that start in 
one of the loops form a infinite tube that ends at the spatial infinity. Now, this tube plays the role of the surface 
of equation ((66)l : the loops Ci and C2 are linked if one of them cuts the tube attached to the other one. 
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V. METRICS g""' Y p^,. IN THE LATTICE 

Now we turn to discrete space. We shall say that a 1-form in the lattice is transverse if its coderivative (divergence) 
vanishes. A projector onto the vector space of l-forms is given by 

55.% = 8,,^6{v -w) + dj^^: (116) 

where (j)™ is a function that satisfies 

9„C = '5(« - w) . (117) 
The following properties ensure that Sl^j^ is a projector as claimed 



^tV- = (118) 
S^^ S^^ = (119) 

5^^ Tf" = T°' \fT/ ST = 0. (120) 



Let us prove (jllSp : 

divSTjw = djvS{v ~w) + djyjdiv(j)™ = -dj^Siv - w) + djw5{v -w)=Q. 
Regarding (|119p . one has 

^T.s^T,-. = E ^Tc.'- s^;^ = E [kcS{v -vn+ 5c." c.] s^;^ 

— "Tjw 2^ 9v" Otjw — '^Tjw ■ 

Finally we prove (|120p : 

The covariant metric in the lattice is a discretization of f[55)) defined as 

g^-J^ = -uj + U- iLi23) ■ (121) 

The following properties of (7™ J™ are " mimetizations" of (|103p . (|102p . and (|100l) respectively 

<?*^^A^ = (7^A)™ (122) 
9™5™^"'=0 (123) 
5'"' =5'^^- (124) 

Let us prove (I122p . which states that the contraction of the metric with a vector yields its curl 

w 

Next we prove that the metric is transverse (property (jl23p ) 

div 5'"^"" = -£ijkdkvdivS{v ~w + L, - ^Li23) = -e^jkdkvdtvSiv ~ w ~ ^Li23) = . 
Also, the metric is symmetric (property (|124l) ) 

gjwiv ^ -Sjif. dkwS{w ~ V +j -iLi23) = -Sijk dkv S{w -V + Lj- ^Li23) 

= -£ijk dkv S{w - w + Lfe + Lj - iLi23) = -£ijk dkv S{w ~ u - + iLi23) = 5™''"' • 
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We take as the discretization of the covariant metric (|94p the expression 

_ ,fc(t>+iLi23-ILj-Lfc) _ ,fc(«+Li-iLi23) fTy^^ 

In order to be a satisfactory covariant metric in the transverse space we ask ^™ to obey the discrete versions of 
equations ([Ml y (P7|. that is 



0™ = -0^"-'^') (without summation in i) (126) 

) = •^^'^^"^ • (127) 

The last equations implies 

d,. K = -5,.,C (128) 

a,.C = 0™ , (129) 

as can be easily demonstrated. As in the continuum, the covariant metric satisfies the properties 

5"'^5m/^ = '^?0 (130) 

ffa/3=5/3a, (131) 

whose demonstration follows. Regarding (|130p one has 

„iv ii \ ' iv av" ^ Pi n c-e-^) J,f>(i'+Li — Lj — Li,) 

H i//jjto — 2-^^ ijav" jw — fciafe<-'fet) i/a(t)+Li-iLi23) jtu ~ t-aikt-ajb 'JkvH'w 

In turn, ()13ip can be established as follows 

_ ,fc(to+Lj-iLi23) _ _ ,fc(t>-Lfc) _ ,fe(t>-Lfc-Lj + iLi23) _ 

yjwiv — fcjifc 'A'l' ^ tjj/c ''^(tu+Lj — iLi23) ~ -''^ ~ yiv jw 

As it happens in the continuum setting, the discrete metrics allow to find the "potential vector" whose curl is 
a given transverse vector field T** 

= (SifA)'' = {4^dAY = (TZA)!' . (132) 

A solution to this equation is given by 

A^^g^.T", (133) 

since 

(TZA)'' = "A, ^ gf" "g, xT^ = Sl^^T^ = T'^ . 

Next we shall study the discretization of the metric gt ^ v associated to the " bundle of paths" prescription defined 
by equation (|114p . Let /i be a continuous path in one of the lattices (say the direct one) starting at the origin and 
ending at the spatial infinity. For each site v of both lattices we define, as before, a path obtained through rigid 
translations of h in such a way that the new starting point be w. It can be seen that the lattice loop coordinate of 
this "bundle of paths" obeys 

;^-(/l^_)=;t'(-+-)(/l^), (134) 

which is the discrete version of ()112p . 

The differential constraint for x™(/ito) takes the form 

d,,x'''{hw)^K'"~w) (135) 

which can be cast as 

S{v -w)= 9„x™(/i-) = E [^"C^-) - x'^''-^'\hj^ = E [^"C^-) - x'"(/»(^+L.))] 
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in virtue of (I134p . From this, we have 

d,^x'''{hra) = ~S{v-w). (136) 
Now we discretize (p'h'y, given by equation piSp in the form 

^ ix™(^») - ^X™(/i.+lJ, (137) 



and define de covariant metric as 



_ ,fc(u+L.i-^Li23) /-ioo\ 

9h ivjw = £ijk Ph w ■ U-JoJ 



^h!w satisfies certain properties that make gt to be a covariant metric metric. They are 

d,,4>h ^a = S{V ~ W) (139) 

(140) 



Vh V — Vh W 



(141) 



To estabhsh (|139p we use the differential constraint as given by (jl35p and ()136p 

= ^9„X™(/1») - \^^vX''"{K+u) = -W)- a„X™(/i.) = 5{V - w) . 

The proofs of (IT40| and plTt use identity (fT34)) 
as weU as 

Hence, we have constructed a covariant metric that mimics, in discrete space, the "bundle of paths" metric of the 
continuum. We write down a final expression for it, that should be compared with (jll4p . and that is obtained by 
substituting ^77)) into (1^5)) 

9H = \e,,u[x^^^+^^-^^^''Hh^) - x'^'^-'^'^HK+u-^uJ) ■ (142) 

VI. SUMMARY AND OUTLOOK 

As we saw, the Gauss Invariant ip{Ci, C2) can be written as 

v(Ci,C2) = x^(Ci)g^,x'^(C2), (143) 

where g^i, is the covariant metric in any prescription. In turn, the Chern-Simons term can be expressed as 

Scs = A^gf^-A, . (144) 

Hence, the covariant and contravariant metrics are the kernels of the Gauss Invariant and of the Chern-Simons term, 
respectively, both in the continuum and in the lattice. This is not the only reason why these metrics are interesting. 
Indeed, they can be seen as building blocks for other knot invariants that are interesting in quantum gravity, as we 
shall briefly discuss. The key object to take into account is the Wilson operator Wa(7), which is the trace of the path 
ordered exponential of the non Abelian Chern-Simons potential 



WAil) = Tr 



P exp j, Aady" 



(145) 
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Witten showed that there is a correspondence between the vacuum expectation value < W{j) > in the non AbeUan 
Chern-Simons theory and the Jones and Homfiy polynomials [IJ . Guadagnini, Martellini and Mintchev [7j found that 
the perturbative expansion of < W{j) > provides a systematic method to obtain analytical expressions of knot and 
link invariants. This perturbative expansion has the form 

/ oo / n ■ \ ri 
dAe''^<''WAh)^N + Nj2[ — ) <W^(7)>^"\ (146) 
n=l \ / 

where the gauge potential A belongs to the SU(N) algebra and k is a coupling constant that serves as perturbative 
parameter. Analytical expressions for the first terms in the expansion are [7|, |9[ 

< W{j) >(i) ^ ^^(^) (147) 

< W{^) >(2) ^ -^l-lp(j) + 1 (< w{^) ' (148) 

< Wi-f) >(3) = - ^(^^~^V (7)+ < W{-f) Wi-r) >(2) -i (< W{-f) >(i))^ , (149) 

where (^(7), p(7) and r(7) are knot invariants. The firs one is the self linking Gauss number, while the other are related 
with the second coefficient of the Alexander-Conway polynomial and the third coefficient of the Jones polynomial 
respectively. In terms of loop coordinates of higher rank [19.], these objects are written as [9|] 

Vil) = X^ffM-x" = 2g,,,,,,x''^^- , (150) 
p(7) = 2{-h,,,,,,,,,x''^>^^''' + 5miM35p.m4X^^^^^^^^) , (151) 

(152) 

where the metrics are taken in the transverse prescription and we have set 

/iMiP2M3 = £"^^3MiaffM2/35p37 and e'^^"^^^ = f dhe''^'6{x - t)5{y - t)5{z - t) . (153) 



These expressions show the important role that the metric plays in the construction of analytical formulae for the 
invariants, as claimed. These invariants are candidates to be quantum states of Loop Quantum Gravity, and it could 
be useful to have mimetic discrete versions of them in order to search for lattice solutions of quantum gravity. A first 
step in this program is the study of the metrics in the space of transverse vector densities that we have presented in 
this article. 
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